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The critical velocity Vc for the onset of quantum turbulence in oscillatory flows of superfluid helium is uni- 
versal and scales as ~ \/ (CCO, where K is the circulation quantum and CO is the oscillation frequency. This 
result can be derived from a general argument based on the "superfluid Reynolds number". Only the numerical 
prefactor may depend somewhat on the geometry of the oscillating object because the flow velocity at the sur- 
face of the object may differ from the velocity amplitude of the body. A more detailed analysis derived from the 
dynamics of the turbulent state gives Ve ~ ^/SKCoJ]5, where j3 --^ 1 depends on the mutual friction parameters. 
This universality is compared with the recently discovered universality of classical oscillatory flows. We also 
discuss the effect of remanent vorticity on the onset of quantum turbulence. Finally, by employing the "super- 
fluid Reynolds number" again, we argue how Vc changes when the steady case CO = is approached. In that case 
Vc scales as k/R, where R is the size of the object. 

PACS numbers: 67.25.dk, 67.25.dg, 47.27.Cn 



I. INTRODUCTION 



The critical velocity for the onset of quantum turbulence 
caused by a macroscopic body oscillating in superfluid he- 
lium is discussed. Experiments with spheres' and vibrating 
wires^ ^ give evidence for a universal scaling property of os- 
cillatory superflows, namely vv ~ y/KCO, independent of the 
particular geometry of the vibrating object that drives the os- 
cillating superflow. This experimental result can be derived 
from the "superfluid Reynolds number" Re^ = vI/k, where / 
is some characteristic length scale"*. Turbulence may be ex- 
pected when Res > 1 . At oscillation amplitudes a = v/co that 
are small compared to the size R of the object, a ^ R, the 
amplitude a is the characteristic length scale while R may be 
considered to be infinite. This leads to the above scaling^. In- 
terestingly, for classical oscillatory flows a universal scaling 
has recently been discovered^ which might be related to our 
case. 

More detailed information can be obtained from the dy- 
namical behaviour of the turbulent vortex tangle as recently 
discussed by Kopnin^ for the case of counterflow turbulence 
at constant velocity. Extending this scenario qualitatively to 
the case of oscillating flows give^ w y^SKCo/ji, where 
j3=A(l — a')— Bais derived from the mutual friction param- 
eters a' and a with coefficients A,Z? ^1 (more precise values 
for A and B would require numerical simulations with some 
particular geometry). Below 1 K where mutual friction is 
small in *He (a', a w 0) we can set j3 = 1 (assuming A, B = 1) 
while towards higher temperatures )3 slowly decreases to 0.79 
at 1 .9 K (the highest temperature where could be resolved 
with spheres). This leads to a slight increase of Vc with tem- 
perature that is in good agreement with the data on spheres^. 
Very close to the Lambda transition j3 might decrease very 
fast and ultimately become negative. This would imply that 
because of large mutual friction superfluid turbulence cannot 
exist anymore, analogous to the situation in superfluid ^He 
above T/Tc^iOM 



The numerical prefactor a/8 2.83 is only an estimate be- 
cause the model is still qualitative as long as there is no rigor- 
ous theory of the dynamical behaviour of the vortex liquid in 
oscillating superflows. Experimentally, it may depend some- 
what on the geometry of the oscillating body, in particular 
when tuning forks are considered, as will be shown below. 
The effect of remanent vorticity on the critical velocity at the 
first transition to turbulence after the measuring cell has been 
filled with helium will be discussed in terms of Kopnin's the- 



ory in Chapter III 



II. UNIVERSALITY 

The dependence of Vc ^ -\/fca) is a universal property of all 
oscillatory superflows and not restricted to ^He alone. To our 
knowledge, however, corresponding data on superfluid ^He 
are not yet conclusive (only 2 frequencies have been studied-^) 
or, in case of Bose-Einstein condensates, are not yet available. 
Very recently, a universal scaling property of classical laminar 
oscillatory flows has been discovered by Ekinci et al.^. These 
authors find that a scaling function exists that depends only on 
the dimensionless product of the oscillation frequency (O and 
the relaxation time T of the liquid. While this scaling applies 
to oscillating classical laminar flows, at present no extension 
to oscillating superflows is available. What we have found- is 
that quantum turbulence occurs only when cox < 1/4 where 
in our case T = 2k/ ^v^ is the relaxation time of the vortex 
liquid in a superfluid flow field Vj as introduced by KopnirPEl. 
Similar to the situation in classical flows^ geometry and di- 
mensions bear no effect. One might ask whether there is any 
relation between these cases. Therefore, it is desirable to have 
a rigorous theory of the dynamical behaviour of oscillating 
superflows. 

An interesting consequence of the frequency dependence 
of vv is that the oscillation amplitude at the critical velocity 
ac = Vc/ CO ^ \J k/g) is equal to the average vortex spacing 
lc = \/ \fLc^ where the vortex line length per unit volume is 
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given h^Lc = {vc/Kf- This is a very plausible result. Fi- 
nally, we note that only the oscillation amplitude determines 
the transition to turbulence and not the classical "Strouhal 
number" cic/R which obviously has no significance for the on- 
set of quantum turbulence^. 



III. REMANENT VORTICITY 

In the experiments both with spherei^^^l and, in particu- 
lar, with wires'^ the critical velocity could be enormously 
exceeded in the first up-sweep of the oscillation amplitude 
and V( could only be determined in the down-sweep after the 
transition to turbulence finally had occurred. We are attribut- 
ing this hysteresis to a lack of remanent vorticity at the be- 
ginning of the experiment. In fact, detailed experiments by 
the Osaka group" demonstrate that remanent vorticity can be 
reduced or even avoided if the measuring cell is filled very 
slowly at low temperatures. In that case no transition to turbu- 
lence was observable up to very high velocity amplitudes of 
1.5 m/s. Starting our analysis with an initial remanent vortex 
line length Lq <C Lc we find from Kopnin 's worlF" that in this 
case the relaxation time is enhanced by a factor (Lf/Lo)'/^. 
(No specific assumptions are made where and how these vor- 
tices are distributed throughout the measuring cell.) From the 
above condition for turbulence we now obtain a larger crit- 
ical velocity v'. ~ O) • /q where /q > lc is the average vortex 
spacing of the remanent vorticity Lq that can be determined 
directly from the measured v[. that is larger than by a fac- 
tor Iq/Ic = {Lc/Loy^^- The critical velocity v'. is reached 
when the oscillation amplitude a'. = v^/o) is comparable to 
Iq: a' ~ lo, which is analogous to the result of Ch.2 where 

^ lc which means that in either case vortex lines must be 
within the reach of the oscillating body for turbulence to de- 
velop. This picture is qualitatively supported by numerical 
simulations^. 




FIG. 1: (Color online) Crossover of the critical velocity Vc{co) from 
oscillatory flow to steady flow (0) at O) --^ k/R^, where R is the size 
of the object. 

For tuning forks the situation is more complicated. Both 
prongs of the fork are oscillating towards and away from each 
other The liquid in between is pushed away or sucked in at 
a speed that will depend on the ratio of the width W of the 
prongs to the spacing D. It is clear that for a ratio W /D > I 
the flow velocity is enhanced by this factor Moreover, be- 
cause of the rectangular cross section of the prongs, the flow 
velocity will be enhanced even further near the edges. Hence, 
the measured values of vv will be lower and will depend on the 
geometry of the forks. Experimentally, there is some evidence 
for that but it is rather qualitative because it is impossible to 
vary the resonance frequency of a fork without changing the 
geometry as welP^. 



V. CROSSOVER TO STEADY FLOWS 



IV. THE NUMERICAL PREFACTOR 

While the experimental remits obtained with spheres, 
wires, and also a vibrating grid'^ follow the above behaviour, 
the data on tuning forks'"* seem to deviate towards values of 
Vc lower than expected. This discrepancy may be attributed 
to a geometrical effect on the prefactor of v^. In general, the 
velocity amplitude of the flowing superfluid varies over the 
surface of the oscillating body. For spheres, e.g., the flow ve- 
locity is largest at the equator, where it is increased by a factor 
of 1.5 compared to the velocity amplitude of the sphere. It 
appears plausible that turbulence originates predominantly in 
this region. For wires the corresponding increase is by a fac- 
tor of 2.0. Experimentally, for spheres a prefactor of 2.8 is 
observed' while for the most detailed data on wires obtained 
by the Osaka group^ the prefactor is 2.1. This is expected 
from the ratio 1.5/2.0 = 0.75 = 2.1/2.8 but the perfect agree- 
ment may be accidental because there is still scatter of the 
data. 



The question arises how the frequency dependence of (ft)) 
connects to the steady case where O) = 0. Of course, a finite 
critical velocity is known to exist in that case too. Clearly, 
there must be a crossover to the steady case before (O goes 
to zero. We again consider the superfluid Reynolds number 
Res = vI/k. At low frequencies the characteristic length scale 
/ must now be the size R of the object because the oscillation 
amplitude a diverges: 7? <C a. This gives then vv(0) ^ k/R, 
which one would also expect on dimensional grounds' and 
which is well known from vortex ring production in steady 
flows (except for a factor l/2;r and a logarithmic correc- 
tion). Therefore, the crossover occurs when (O ^ k/R^ (see 
Fig.l), i.e., when the average vortex separation lc at the criti- 
cal velocity or, equivalently, the critical oscillation amplitude 
Oc ~ lc ^ \J k/co begin to exceed R. For our spheres of size 
(diameter) R = 0.2 mm we obtain Vc(0) ~ 0.5 mm/s and a 
crossover at o ~ 2.5 s^' which is much smaller than the os- 
cillation frequencies of our spheres > 750 s^' and, therefore, 
could not be observed in our experiments. For wires, how- 
ever, having a size of only 3 micrometer^ the transition will 
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occur already at 1.1 • 10"* s^^, or 1.8 kHz. For frequencies that 
are substantially lower, the critical velocities will reach the 
constant level of Vc(0) ~ k/R = 33 mm/s (when numerical 
prefactors are taken into account these numbers will change 
by a factor of order 1). 
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